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This article is the first one in a series of two dealing with the concept of "resonant metalens" we recently 
introduced [Phys. Rev. Lett. 104, 203901 (2010)]. Here, we focus on the physics of a medium with finite 
dimensions consisting on a square lattice of parallel conducting wires arranged on a sub-wavelength scale. This 
medium supports electromagnetic fields that vary much faster than the operating wavelength. We show that such 
modes are dispersive due to the finiteness of the medium. Their dispersion relation is established in a simple 
way, a link with designer plasmons is made, and the canalization phenomenon is reinterpreted at the light of our 
model. We explain how to take advantage of this dispersion in order to code sub-wavelength wave fields in time. 
Finally, we show that the resonant nature of the medium ensures an efficient coupling of these modes with free 
space propagating waves and, thanks to the Purcell effect, with a source placed in the near field of the medium. 

PACS numbers: 41.20.-q, 81.05.Xj, 78.67.Pt 



I. INTRODUCTION 

In a recent papei^, we have introduced the concept of "reso- 
nant metalens" which can be used to break the diffraction bar- 
rier. We defined it as a cluster of coupled resonators, arranged 
on a subwavelength scale. When illuminated with broadband 
electromagnetic waves this macro-resonator radiates waves 
that contain information of an object placed in its near field. 
By measuring this radiation one can recover sub-wavelength 
details of the object, and thus focus (or image) waves on a 
deep subwavelength scale, not limited by the Rayleigh crite- 
rion. 

In the original papei^, we used metallic wires, half a wave- 
length long, as sub-wavelength resonators and we propose to 
explain the physics of this example. We focus on the reso- 
nant metalens introduced initially ^ consisting of a square lat- 
tice of N X N (N = 20) identical metallic wires, aligned 
along the vertical axis (defining the longitudinal z-direction). 
The length of the wires along z is equal to 40 cm and their 
diameter to 3 mm. The period lattice in both transverse direc- 
tions (xy plane) is 1.2 cm, which is roughly equals to A/ 70 
for the first intrinsic resonance of a single wire. 

The lens of interest here falls within the class of "wire 
medium". We define the latter as a uniaxial medium formed 
by a periodic lattice of conducting wires with small radii com- 
pared to the lattice periods and the wavelength. In microwave 
applications this medium is often seen as an interesting arti- 
ficial dielectric, which presents the great advantage to have a 
negative effective permittivt>^. Some attention to wire me- 
dia has been paid in the realization of left-handed media as 
composite media made from lattices of long conducting wires 
and split ring resonators^ . Usually, quasistatic models for this 
medium are carried out for wave propagation perpendicular 
to the wires which show that, for electric-field polarization 
along the wires, the medium is characterized by a frequency- 
dependent effective dielectric constant. 

Belov and his colleagueP^lhas studied this kind of medium 
in details, and especially the propgation along the wires. They 
showed that this medium can be described by means of an 
anisotropic dispersive dielectric tensor. Also, they demon- 
strated that the wire medium supports propagating modes, so- 



called transmission line modes, which travel along the wires. 
During the whole article, we will focus only on these modes 
which explain our initial result^ and permit to give a simple 
explanation for our meaning. 

These modes which propagate in an adispersionless way in 
an infinite wire medium gives rise to dispersion when adding 
finiteness along the wires. While in the original papei^ we 
insisted on the fact that the dispersion is a key issue for sub- 
wavelength focusing/imaging with our method, we will give a 
simplified model for the physical mechanism responsible for 
it and a discussion justifying our simplification is presented 
in appendix. We show why such a medium can be used to 
code a sub- wavelength wavefield into a temporal/frequencial 
signature thanks to its dispersion relation. 

This will lead us to note the link between our work and 
Pendry's designers plasmons^, hence introducing the waves 
propagating in a finite length wire medium as new spoof plas- 
mons. We will as well comment the very interesting effect 
which occurs at a specific frequency, known as canalization 
regimJ^, at the light of our simplified approach. 

The other key issue for a resonant metalens is the far field 
radiation of the sub- wavelength features. Usually, a sub- 
wavelength wave field is considered as evanescent, which 
forbids its propagation toward the far field. But, the notion 
of evanescent waves, which is responsible for the diffraction 
limit, is a mathematical formalism which only fits that of in- 
finite interfaces. Using finite transverse dimensions and end- 
ing on the medium presented initially ^ , we show that the sub- 
diffraction details of the latter actually contributes to the far- 
field radiation. 

Eventually, it is commonly admitted that their contribution 
to the total radiation is much weaker than those of propagat- 
ing diffraction limited waves. Due to the resonant nature of 
the medium of interest, another physical mechanism which 
is well known from opticians intervene: the Purcell effecf^. 
This effect increases the coupling of the source with radiating 
modes, and surprisingly it counterbalances the effect of the 
weak radiation of the sub- wavelength features. 
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II. THE INFINITE WIRE MEDIUM: ADISPERSIONLESS 
PROPAGATION OF TEM BLOCH WAVES 

We consider an infinite "wire medium" consisting of infi- 
nite Perfect Electric Conductors (PEC) wires afigned along 
the 2;-axis, forming in the xy-plsmc a two dimensional square 
lattice (Fig. [T]). We assume a thin wire approximation, ie. 
their radii are really small compared to the lattice constant. 
We are looking for electromagnetic fields {S;l-L) solutions of 
the following Maxwell problem and satisfying the boundary 
conditions: 
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V xn = eo 
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where fiQ, Sq are respectively the vacuum permeability and 
permittivity. Here we have assumed that the electromagnetic 
fields will propagate in vacuum only: the presence of metal- 
lic surfaces will be taken into account through boundary con- 
ditions only. Using the invariance of the medium along the 
z-axis and choosing a time dependence in e~*^^ we define 
time-harmonic two-dimensional electric and magnetic fields 
E and H as: 



n{^±,z,t) = J?e (H(x^)e^(^^-^')) (2) 

where stands for the coordinates in the {xy) plane and 
k is the propagation constant of waves along the z-axis. E 
and H are complex valued fields depending on two variables 
(coordinates x and y) but still having three components along 
the three axes. The medium of interest is an uniaxial medium 
and we propose to decompose the electromagnetic fields with 
respect to the geometry of the medium: 



E(x^) =E^(x^)+^ll(x^)e, 
H(xx)=H^(x^)+i7||(x^)e, 



(3) 



Eventually, due to the periodic nature of the medium in 
the transverse plane, the problem becomes the well known 
problem in solid state physics which consists in looking for 
Bloch waves solutions /7k^ (where U stands for E^, H^, E\\ 
or H\\) that have the form: 



[/k^(x^) = e 



ik_L .X_L 



(4) 



where is the so-called Bloch wave number (or quasi- 
momentum in solid state physics) and t/(xx) is a function 
that satisfies the periodicity of the wire medium (U{'k_\_) = 
U (xx + R) where R is a lattice vector). With these geomet- 
ric considerations the Maxwell system reduces to: 
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FIG. 1. (top) The geometry of the infinite wire medium: an in- 
finite square array of parallel perfectly conducting thin wires. We 
focus on the frequency range where the lattice parameter is sub- 
wavelength, (bottom) Dispersion curves for the TE and TEM types 
of Bloch waves supported by the medium. The TEM Bloch modes 
are always propagating with a flat dispersion curve. For the geometry 
under consideration (r <^ a <^ X) the TM Bloch modes are always 
evanescent^. 



iuj/ioH\\ez = (^i^± + x Ex 
-zcjeoEx = ez X (^ikil± - (^k± + V ) i^n ) (5) 
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Then, we study independently the three types of Bloch 
modes that can exist with respect to the wires: the transverse 
magnetic (TM) where H\\ = 0, the transverse electric (TE) 
where = 0, and the transverse electromagnetic (TEM) 
modes where H\\ = =0. The final equation to be solved 
for each type of mode becomes: 



TM 
TE 
TEM : 



k^- 



■(^)')^l|-2z(kx.V)^l|-A^II=0 
{^f) ^||-2z(kx.V)^||-A^||=0 
(fe2_(c^)2)Ex=0 (6) 



Solving the boundary value problem for TE and TEM 
modes is easy since they do not have any electric field compo- 
nent along the wires, and a dispersion law is easily obtainable 
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from the previous equations. For the TM modes this boundary 
problem is a tedious work which has already been performed"^. 
As a summary, the solutions of the Maxwell eigenproblem in 
the first Brillouin zone must satisfy the following dispersion 
relations (Fig.[T]): 

TE ■.k^^ + e = '^ 

TEM :k^ = % (7) 

where kp is the so-called plasma wave number, a parameter 
that depends on the lattice parameter a and the wires radii r 
and its expression can be found in the litterature'^. 

From now, we focus on the range of angular frequencies 
uj where the associated vacuum wavelength is really higher 
(> 10) than the distance between two wires. In that case 
and in the thin wire approximation, the plasma wavenumber is 
greatly highei^ ^ than the free space one and the TM waves are 
always evanescent. The propagating TE modes along the z- 
direction are included inside the vacuum light cone. And the 
TEM Bloch modes are always propagating with a flat disper- 
sion curve independent of the norm of the Bloch wave number 
k±_ . Clearly, for the large k±_ corresponding to sub- wavelength 
features, the TE and TM modes are evanescent, while the 
TEM ones are propagating. We will only consider the propa- 
gating modes for this range, ie. the TEM Bloch modes. Thus, 
the fields inside the structure for a given k±_ can be described 
by a unique scalar potential (j){x±_)'. 



III. THE FINITE WIRE MEDIUM ALONG THE WIRES 
AXIS: A DISPERSIVE MEDIUM 

A. Snell's law for a wire medium-air interface 

Unfortunately, the infinitely extended wire medium cannot 
exist and introducing a boundary condition leads to interesting 
results. Here, we consider a TEM Bloch wave with transverse 
wave number = /c^ex traveling in the half space z < L/2 
filled out with a wire medium. This Bloch wave encounters 
a boundary with vacuum at z = L/2. The incident angle 
6i is defined by k^/k. This wave generates a reflected wave 
with angle Or and a transmitted wave with angle Ot. The first 
step consists in finding an equivalent of the Snell's law for 
this interface. Considering the conservation of the tangential 
component of the wave number at the interface (here k^) and 
introducing the dispersion relation on both sides of the equal- 
ity we obtain: 

Of = — 9i 

tan6>i = sin6>t (9) 

From this refraction's law the critical incident angle Oc = 7r/4 
appears: when the incident angle is higher than 9c (ie. for 
evanescent waves in the air region) a total internal reflection 
is observable. 

In order to simplify the calculation for continuity of tan- 
gential components, as explained in the previous part, we will 
neglect the evanescent TM wave that is generated in the wire 
medium at the reflection and the high diffraction orders. A 
discussion concerning those approximations is presented in 
appendix. Then, according to equation ([S]), the electric fields 
on both sides of the interface write: 



Ek^(x^) = -W^^-"^^(x^) 

. _ ^ iik±x-\-k(z-L/2)) 

Hk^(x^) = X W^^-"^^(x^) (8) 



E^^) = ^^Qe^(fc^^-fc(^-^/2)) (10) 




^(xx) must satisfy the periodicity of the medium as well as 
the boundary conditions at metallic surfaces. This potential 
can be seen as a spatial envelope that permits to verify the 

constraints, and it multiplies a plane wave, as a consequence E*^^^ = tE^e^^^^^^'^^^'^^'^^^ 

of the Bloch theorem. 

During the rest of this article, we will always neglect the 
high diffraction orders due to this potential and consider it as ^nd the associated magnetic fields write: 
a constant potential. We will always talk about plane waves 
instead of Bloch waves neglecting the Bloch potential. One 

important consequence of this hypothesis is that we can arbi- _ ^i{k±x+k{z-L/2)) 

trarily change the coordinate axes without lack of generality, Moc 
and especially we will often choose — ^_l6x. 

Another interesting result concerning the TEM Bloch tEq i(^k±x-k{z-L/2)) 

modes is that the group velocity Vg = Vk(^) = cGz is ~ jj^qc^ 

constant and oriented along the z-axis. Therefore any spatial 
wavepacket in a transverse plane propagates without distor- 

tion along the 2:-axis and one can use this property for guiding H*^^^ = — 9.^<k±x+Kiz-L/2)) [ i | (H) 

and imaging sub wavelength features 
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where k = ^ stands for the free space wave number, and k = 
— k±'^ if Oi < Oc or = i\^k±'^ — k'^ if the total internal 
reflection occurs. Applying the continuity of the tangential 
components of the fields at the interface and introducing the 
refraction law of equation ([9| yields to the following reflection 
coefficient: 



(12) 



Interestingly, in the peculiar case of total internal reflection, 
r becomes a complex coefficient with unity norm, and the re- 
flection induces a phase shift: 



(f{k±^ k) = 7T — 2 arctan 



(13) 



This phase shift is a well-known phenomenon for total internal 
reflection and one of its consequences is the Goos-Hanchen 
shifP along the x-direction. Here, this phase shift has an 
importance when adding a second boundary at z = —L/2. 
Now, the wire medium with finite dimension consists of wires 
aligned along the z-direction with the same length L (Fig. [2]). 
With this new constraint the field inside the medium experi- 
ences a Fabry-Perot like resonance, or equivalently a guided 
propagation along x-axis, when the optical path length S for a 
round trip satisfies: 



2kL^2ip{k±,k) = 2n7r 



(14) 



where n is an integer. With equation ([13]) and ([14]) we can 
write a dispersion relation for the guided waves inside the 
bounded wire medium for the first resonance (n=l)\ 



tan k 



Vk^^ - k^ 



(15) 



This implicit theoretical dispersion relation is plotted in 
terms of the first Fabry-Perot frequency fm versus kj_ (both 
normalized to /o and ko corresponding to an ideal Fabry-Perot 
of length A/2) in Figure [2] The asymptotes of the light line at 
low frequencies and the intrinsic resonance frequency at large 
values of k± are shown. 

Rigorously, the TM evanescent waves as well as the high 
diffraction orders (and thus the TE diffracted waves inside the 
wire medium) induced at the reflection adds small corrections 
to the phase shift. This would lead to a more complicated dis- 
persion law that depends on the wire medium parameters. For 
simplicity and because this dispersion is in a sufficiently good 
agreement with simulation results (Fig. [7|, we have preferred 
not to introduce those waves. More sophisticated dispersion 
laws are presented in the appendix. We believe our calcula- 
tion that only takes into account the evanescent field at reflec- 
tion is the easiest way to understand why the z-bounded wire 
medium introduces the dispersion. Because this dispersion is 
a key issue for the concept of resonant metalens, we have cho- 
sen to make those approximations for the scope of this article, 
but obviously one can add small corrections to this law. 



B. Fields profiles inside the medium 

In order to describe precisely these propagating waves in- 
side the wire medium we consider that the region between 
z = ±L/2 is filled out with a wire medium (Fig. [2]). We 
know that the field inside this region is TEM with respect to 
the wire but the invariance along z is broken. Still considering 
the first order of diffraction, neglecting TM evanescent waves, 
and choosing the x-axis in order to have = k±ey^, we can 
express the fields inside the wire medium as: 



ETEM = e^^-^- ( 








Htem = ^e^^^-^^ dEo/dz 



(16) 



And, from the propagation equation ([5]) of the TEM Bloch 
modes inside the wire medium it follows: 



Eq{z) = Ci sin (^^^ + C2 cos 



/UJZ 

V c 



(17) 



In the air region on both sides of the wire medium the elec- 
tromagnetic fields decompose onto a plane wave that we can 
write using the previous notations: 



±n/k 

-k^/k 



El 







(18) 



Then applying the continuity of the tangential fields at the 
interface we reach the matrix form equation: 



/ sin(/cL/2) 
-cos(/cL/2) 
- sin {kL/2) 

\- cos {kL/2) 



cos {kL/2) 
sin {kL/2) 
cos {kL/2) 
-sm{kL/2) 
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(19) 



The system has a non-zero solution only if its matrix de- 
terminant is null. The first solution of this cancellation gives 
another expression of the dispersion relation of equation ([HI: 



cos 



_k_ 



(20) 



This expression is actually the same as the one presented 
previously in equation ( [T5] ), and one can develop the sine func- 
tion from its cosine and recover the tangent function. From 
this expression, we can eventually express all of the electro- 
magnetic components. In the wire medium region it gives: 
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FIG. 2. (top) The geometry of the infinite wire medium with finite 
length along the wires, (bottom) The dispersion relation of the first 
resonance frequencies fm (normalized to the original Fabry-Perot 
frequency /o) versus the transverse Bloch wave number k± (nor- 
malized to the original Fabry-Perot wave number ko) from equation 
|T5|l,([20| or ( [25| l. The asymptotes of the light cone and the original 
resonance frequency are shown (dashed-lines). The inset shows how 
the three ways of calculation (reflection coefficient, guided waves 
expression and effective Fabry-Perot length) leads to the same result. 
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C. Energy density considerations 

However, we can use a more intuitive approach to deter- 
mine the dispersion law. In fact, depending on the trans- 
verse wavevector k±, the Bloch waves penetrate in air up to 
a distance whose inverse is the quantity previously introduced 

VkA_^ - k^. Again, a more sophisticated model would add 
correction to this penetration depth taking into account high 
diffraction orders and TM waves (see appendix). Each propa- 
gating mode in the structure experiences the Fabry-Perot res- 
onance, and the length of the cavity is not simply the length L 
of the wire medium, because of the penetration of the modes 
inside air. A common way to describe the effective length of 
the real cavity is to calculate the quantity: 



u{x^ z)dz 



/CO 
-oo 



(23) 



z)rdz 



where u{x , y ^ z) is the energy density of a given mode. Us- 
ing the fact that the modes are TEM inside the structure, the 
density of energy inside the structure is constant along the 
2:-direction, from —L/2 to L/2. For the sake of simplic- 
ity, we set this value to 1, which does not hamper our de- 
scription since this value simplifies in Leg. The energy den- 
sity decreases exponentially at the interfaces z = —L/2 and 
z = L/2 and taking into account the three non-zero compo- 
nents of the fields, the calculus simply gives: 



L. 



(24) 



This expression has merit to show explicitly the increase 
of the Fabry-Perot cavity length due to penetration depths at 
interfaces. We finally consider that the resonance condition is 
obtained if Lgff = nX/2 for the nth Fabry-Perot like condition, 
which leads after rearrangement to the dispersion relation ini- 
tially presented^ for the first Fabry-Perot resonance: 



One can notice that for the first resonance the electric field 
is maximum at both ends of the wire medium, while the mag- 
netic field is maximum at half the wires length, and is near 
zero at both ends (fig. [2]). And in the vacuum regions for high 
Bloch wave numbers (ie. the waves are evanescent in air) it 
gives: 



fo 

f 



ko 



(25) 



where ko and fo stands for the original first Fabry-Perot Uke 
resonance, namely ko = n/L and fo = c/2L. 



Hoc k±_ 



(22) 



D. Link with designer plasmons 

The dispersion curve deserves some comments. Indeed, 
this relation is similar to the ones characterizing surface plas- 
mons polaritons (SPPs) and we have demonstrated the exis- 
tence of guided modes inside the wire medium (with finite 
length) with high wavenumbers. But, at low frequencies, 
where the limit of perfect conductor is applicable, it is well 
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known that perfect metals do not support electromagnetic sur- 
face modes, forbidding the existence of SPPs. 

However, it has been demonstrated theoreticall}0 and 
experimentally that bound electromagnetic surface waves 
mimicking SPPs can be sustained even by a perfect conductor, 
provided that its surface is periodically corrugated. Two dis- 
tinct geometries of modulation of the flat perfectly conduct- 
ing surface have been introduced^^: a one-dimensional array 
of grooves and a two-dimensional hole array. The frequen- 
cies of the supported modes scale with the geometrical size 
of the corrugations in the perfect conductor approximation. 
Maier and his colleagues also demonstrated the existence of 
pseudo surface waves in the low frequency regime: structures 
consisting of closely spaced metal rods in one dimension can 
guide electromagnetic waves. In the present case, the same 
configuration has been investigated but in the two dimensional 
case, and we have demonstrated the same guiding properties. 



veloped here. From previous calculations, introducing equa- 
tion (15) into equation ([12]) leads to the following reflection 
coefficient: 



Interestingly, the equation ( 15 ) is actually the same disper- 
sion law as the one obtained for the spoof plasmons^^. Again, 
in both cases the high diffraction orders and the mode conver- 
sions are neglected when calculating the dispersion law, and 
this is the reason why the two laws are the same. Because of 
the matching of the dispersion laws at the first order, one can 
imagine to couple the two kinds of structure. Similarly, the 
thickness of the layer that supports these plasmon-like waves 
scales with the operating frequencies since the height of the 
wires is near A/2. But, here, only the length of the wires has 
an influence on the resonance frequency. It is not the case for 
the array of grooves where the groove's width and the distance 
between two of them have an impact on the dispersion curve. 
In the array of pinholes, the fundamental mode inside a hole 
is considered which implies a dependence of the hole dimen- 
sions in the dispersion. Here, the distance between wires as 
well as wires diameters do not influence the dispersion law 
and the guiding properties of the medium do not seem to be 
significantly altered when wires are displaced from their orig- 
inal positions. Another interesting property is that the wire 
medium is essentially made out of vacuum, thus losses in- 
duced by metallic surfaces are less important than in the cor- 
rugated surfaces or in the surface with pinholes. 



E. Link with canalization regime 

Eventually, we would like to add comments concerning the 
results obtained by Belov and his colleagues^. They have 
used the wire medium for performing a pixel to pixel imag- 
ing with a resolution thinner than the diffraction limit. The 
same near field lens has also recently reached the near infra- 
red spectrum^^ and sub-wavelength imaging has been proved 
experimentally using gold nano- wires. In the two cases, the 
approach is monochromatic and the operating frequency of 
their lens corresponds to the intrinsic resonance of a single 
wire, namely /q. At this exact frequency the field distribution 
at an interface is recovered at the other interface, and by scan- 
ning near field of one interface they perform sub- wavelength 
imaging. We would like to prove that their results correspond- 
ing to the canalization regime are consistent with the ones de- 



1 — i tan(7r/L/c) 
1 tan(7r/L/c) 



(26) 



Therefore, at the intrinsic resonant frequency of a single 
wire /o this reflection coefficient is equal to — 1 and does not 
depend on the Bloch wavenumber k± since the dispersion re- 
lation has no solution at this frequency. Because this coeffi- 
cient is real, the reflection does not imply any phase shift and 
especially the Goos-Hanchen shift is null. As we have seen 
before, the group velocity for the TEM Bloch waves inside the 
wire medium is oriented along the wire axis, so if there is no 
Goos-Hanchen shift at the reflection, there is no propagation 
in the transverse plane. At /o, whatever the Bloch wavenum- 
ber k± considered the reflected wave propagates along the 
wires as well as the incident one, hence giving rise to the so 
called canalization regime. This regime implies that the elec- 
tric field has the same profile at both interfaces. When placing 
a sub-wavelength source at an interface they managed to im- 
age it half a wavelength away at the other interface with a 
resolution far better than the diffraction limit^. 



F. Coding of sub- wavelength details 

The dispersion is a key issue for the concept of resonant 
metalens we introduced^ The importance of this dispersion 
manifests when a sub-wavelength source is placed near one 
of the two interfaces. Actually, the two interfaces correspond 
to the object plane of the lens and the lens permits to obtain 
the field profile inside one of these transverse planes. This 
profile decomposes with a unique set of amplitude and phase 
on the modes supported by the wire medium. Typically, it can 
be expressed in terms on k± since it corresponds to a Fourier 
transform: 



E{^±) = JJ £:(kx)e^^^-^^dkx (27) 

where in the present case the integration domain is limited to 
the first Brillouin zone of the lattice, that is to say when one 
of the component of the Bloch wavenumber reaches ir/a. 

This formalism makes sense when considering broadband 
illuminations. Interestingly, from previous calculations and 
ignoring the degeneracy, a given Bloch wavenumber k± is 
associated to its own angular frequency uj. For the vacuum 
dispersion relations this is also the case but here k± can be as- 
sociated to a spatial scale that is deeply sub- wavelength com- 
pared to the associated vacuum wavelength. This property 
permits to transform the function £^(k^) into a function that 
depends only on frequency: this is the frequency coding of 
the sub- wavelength details. This property means that, when 
the medium is illuminated from one of the two object planes, 
registering the frequency spectrum of a wave packet gives the 
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spatial Fourier transform of the source. The spatial features of 
the source can now be expressed in terms of uj: 

^(kx) = H^) (28) 

Eventually, it is well known that time and frequency are two 
distinct representations of the same phenomena. A broadband 
illumination means the emission of a short pulse, and the fre- 
quency signature J^{uj) for the sub- wavelength features means 
a temporal one f{t). The imaging procedure associated to this 
concept is the scope of the next article, but at this time one can 
easily understand that the resolution must be limited by the 
distance between two wires. By exploiting the dispersion of 
resonant modes with broadband illuminations we can record 
the sub- wavelength details of an object. In the next part we 
will see that this information can be carried toward far field. 




FIG. 3. The geometry of the wire medium with finite dimensions in 
the 3 spatial directions. Due to the finiteness this medium radiates 
field toward the far field and superimposed on the medium we have 
represented a sketch of a radiation pattern for the first resonance. 
This is just an example since the next figure shows that the directivity 
pattern in the {xy) plane can be different. 



IV. A WIRE MEDIUM WITH FINITE TRANSVERSE 
DIMENSIONS: FAR FIELD RADIATION OF 
SUB-WAVELENGTH FEATURES 

Until now, everything happened inside or in the near field 
of the wire medium. Now, we consider the previous wire 
medium but with finite dimensions in both transverse direc- 
tions (Fig. [3]). Introducing this finiteness has a first conse- 
quence which is the quantization of the transverse wavenum- 
bers. In solid state physics, the quantization is done by intro- 
ducing Born Von Karman boundary conditions'^ which con- 
sists in cycling a crystal. When studying crystals the number 
of atoms considered is large enough to allow this approxima- 
tion. In the present case, the number of wires in both trans- 
verse directions is only 20 and this boundary condition is not 
well suited. 

But, introducing finiteness results in the same quantization 
effect. At each interface, the plasmon-like wave traveling in- 
side the wire medium encounters a boundary with vacuum 
and gives birth to a reflected one. The analytical expression 
of the reflection coefficient is not easy to calculate since the 
field in the air region is the superposition of many plane waves 
due to the finiteness along the wires. Then, the interferences 
inside the wire medium are constructive only for a discrete 
spectrum of Bloch wavenumbers which characterize the sta- 
tionary eigenmodes of the system. Thus, due to the finite- 
ness of the array of wires, the transverse Bloch wave num- 
bers are quantified: = §(m.ex + n.Gy), with integers 
(m, n) G |0; A/' — Ip and D the size of the medium in the 
X and y directions, D = a{N — 1). Again, from the previ- 
ous section, it is important to keep in mind that each /cx(m,n) 
is associated to its own resonant frequency /(m,n)' ensuring a 
frequency (or equivalently a time) coding of spatial features. 

The second consequence of these interfaces is a conver- 
sion from the plasmon-like modes to propagating waves in 
the free space region. Intuitively, with an approach of coupled 
resonators, one easily understands that each wire generates a 
monopolar radiation (for the first resonance) which is essen- 
tially z-polarized. The wire medium can therefore give rise to 
a superposition of monopolar radiations. With the approach 



we used until now, the simplest way to explain the radiation 
comes from the fields profiles inside the wire medium (inset 
of figure [2]). We have already seen that the fields cannot exit 
the structure from the top or the bottom interfaces because 
the waves are evanescent. The only way to escape the wire 
medium is from lateral interfaces. Since the x and y com- 
ponents of the electric field present an odd profile for the first 
Fabry-Perot resonance, they cannot contribute to far field radi- 
ation. While the z-component of the electric field can be mod- 
eled by two small dipoles pi^2 positioned at z = ±L/2. The 
field radiated by the superposition of these two infinitesimal 
dipoles writes, in the dipolar approximation and in spherical 
coordinates: 

^ ^ M0^% sin <l> ^^fe^ / i^sin<^> ^ ^i(-^sin<^>+(^)\ 

47rr V J 

(29) 

where ^ corresponds to the zenith angle, and (p stands for 
the phase difference between the two dipoles. Depending on 
the parity of the Fabry-Perot resonance order the two dipoles 
are in phase or out of phase {ip = or tt). Typically, for the 
first Fabry-Perot resonance, or for the field profile presented 
in figure [2j the two dipoles are in phase, ie. (p = 0. Thus, 
for the first Fabry-Perot resonance, the radiated amplitude is 
maximum in the plane z = and it can be studied as a two 
dimensional problem. 

Now we consider that we have such small dipoles on both 
interfaces of the wire medium with a 2D spatial distribution 
given by the z-component of the TEM Bloch mode at the in- 
terfaces. We perform the projection of this mode onto the 
basis of 2D free space plane waves. The field at the inter- 
faces of the structure has the periodicity given by the Bloch 
wavenumber and has the form: 



E^^ (xsm\^k^^{x - — )j sm\^k^^{y- — )j Uoix.y) 

(3( 

where Hi) is the 2D rectangle function of width D. 
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FIG. 4. The two dimensional projection of three distinct Bloch 
eigenmodes onto the free space waves. The left colormaps show 
that the finiteness of the wire medium in the transverse plane lead to 
a superposition of 4 sine functions. The right ones is the projection 
of the left ones onto the light cone: when considering object with 
finite dimensions a certain amount of energy can reach the far field. 
These colormaps also show that this representation gives information 
about the directivity pattern of an eigenmode. Depending on its k± 
monopolar, dipolar or quadrupolar radiation can occur. 



In the Fourier domain, this expression transforms into a 
sum of four sine functions positioned at the corners of a rect- 
angle of width 2k±^ and length 2k_\_y (Fig. |4]). Such a repre- 
sentation contains all information about the propagating fields 
generated by the sub- wavelength modes. The intersection of 
the latter with a disk of radius uj/c (green line in figure |4]) 
gives the efficiency of the radiation as well as the directivity 
pattern. 

Averaging this intersection along the disk resides in an es- 
timation of the radiated energy in the far field for a given 
mode. The full calculation of the radiated energy resides 
in a projection onto a set of Hankel functions and is quite 
hard to perform, but an estimation of this quantity can be 
done. Due to the fact that we sum 4 sine functions whose 
operands are proportional to Dk±^y, and due to the fact 



that the transverse dimension is lower than the operating 
wavelength, the radiation's efficiency in energy is roughly 
equal to l/{167r).{DLk±^k±y)'^: the smaller the transverse 
wavenumber the higher the radiated amplitude. Fortunately, 
we will see later that the Purcell effect compensates for this 
low radiation efficiency. 

In the same time, the spectral representation of the modes 
gives information about the directivity pattern of the radiated 
field. The intersection between this spectrum and the circle 
of radius k directly gives the radiation pattern of the modes. 
A point on the circle with coordinates (cos 6>,sin 6>) permits 
to obtain the radiated amplitude of the mode in the direction 
defined by the angle 0. The directivity pattern is therefore 
fully described by this representation. It is important to keep 
in mind that this representation is in the (x?/) -plane and that 
the radiation in a plane containing (Oz) is like the one of a 
single dipole antenna. To underline this aspect we represent 
the projection of the previous /c- space spectra onto this circle 
(Fig.g. 

Depending on the mode parity, or equivalently of the parity 
of m and n, the Fourier spectrum has to keep the same parity. 
For a given direction in the transverse plane, if the integer 
describing the wave number component is odd, the spectrum 
is even: it results in a null a value for the spectrum at the 
origin. On the contrary, an even integer does not impose a 
null value at the origin, and it does not cancel over because in 
the present case the transverse dimension D is lower than half 
a wavelength. Taking into account the two directions it results 
in four directivity pattern possibilities: 

• the odd-odd modes give monopolar radiation; 

• the odd-even (resp. even-odd) modes give a dipolar 
along X (resp. y) radiation; 

• the even-even modes give a quadrupolar radiation. 

Therefore the resonant metalens possesses four spatial de- 
grees of freedom, representing as many information channels 
exploitable for imaging and focusing At a frequency where 
some modes are degenerate, registering far field in different 
directions allows a discrimination between modes, by the way 
increasing the number of subwavelength information regis- 
tered. 

This type of directivity patte rns can also be interpreted 
in terms of corner radiatiorP^*^ which have been introduced 
in acoustics. In calculating radiation from vibrating square 
plates, a specific kind of modes, sometimes termed corner 
modes for the reason that only the corner quarter-cells con- 
tribute significantly to the radiation, have been described. 
Those modes satisfy the condition k < ttitt/D and k < 
nn/D. These modes correspond to subsonic modes inside 
plates, analogue to the Bloch eigenmodes considered here. 
If we decompose the field onto subwavelength secondary 
sources, most of those sources tend to cancel one another. By 
grouping secondary sources onto dipoles or quadrupoles (fig- 
ure [5]), it is shown that the effective radiating area is limited 
to the four corners of the square, and the type of directivity 
pattern is easily obtained. 
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FIG. 5. Illustration of corner modes radiation, (left) The z- 
component of the electric field profile for a given Bloch eigenmode 
at one of the two interfaces z = ±L/2. (right) Decomposition onto 
secondary sources demonstrating that the radiation is due to the cor- 
ners. In this case one easily understand that a dipolar radiation is 
obtained. 



V. PURCELL EFFECT 

The idea we introduced in the original papei^ was to use 
the radiation of these subwavelength features to image an ob- 
ject placed near the resonant medium. Thus the last issue to 
be resolved is the coupling between the object and the lens 
itself. Spontaneous emission is not only an intrinsic property 
of the atoms or molecules in a specific material but is also 
governed by the electromagnetic properties of the host. In the 
1950s, E.M. Purcell^^ demonstrated that the rate of sponta- 
neous emission for a light source depends on its environment: 
it can be increased when placing the source in a resonant cav- 
ity. This phenomenon is nowadays well known as "Purcell 
effect". 

Later, E. YablonovitclPD predicted theoretically that a "pho- 
tonic" material can control the rate of radiative recombination 
of an embedded light source. This led to the huge domain of 
photonic bandgap materials where this effect is used for in- 
hibiting the spontaneous emission of light. 

Also, many theoretical works concern the effect of an ab- 
sorbing dielectric on spontaneous emission and level shifts of 
an embedded atom, using a Green's function approach^^ or 
using a macroscopic Lorentz's cavity model^^, or using equa- 
tions of motion for a collection of two-level atoms embedded 
in a dielectric host medium^^. As an atom in a dielectric is 
surrounded by many other two level resonators, these works 
are similar to the concept of resonant metalens we study here. 
Interestingly, we have already discussed the frequency shift 
in the case of the resonant metalens and here we discuss the 
increase of the spontaneous emission rate. Namely, in the res- 
onant metalens case, the Purcell effect implies an increase of 
the coupling between the source and the wire medium since 
it acts as a cavity surrounding the source and modifies the so- 
called local density of states (LDOS). 

The computation of the eigenmodes readily yields the qual- 
ity factor of the modes resonance, a factor influencing the 
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FIG. 6. Measured return loss for a small dipole with (blue) and with- 
out (red) the wire medium. Due to the resonant nature of the wire 
medium a completely reactive dipole becomes radiative, and this ef- 
ficiency grows with the resonance Q factor, and thus with frequency 
here. The discrete nature of the blue curve still demonstrates that the 
subwavelength features are coded in the field spectrum. 



spontaneous emission. From the previous calculation, and 
in the deeply sub-wavelength limit, we showed that the ra- 
diation's efficiency is roughly proportional to {DLk±^k±y)'^: 
the lower the transverse wave number, the higher its radiated 
energy. In other words the resonance quality factor, which is 
the ratio between the energy stored in the resonant medium 
and the radiated energy, grows as {DLk±^k±y)'^ . As stated 
by E.M. Purcell it results in a better rate of spontaneous 
emissions for the highest transverse wave numbers. 

This effect manifests itself by a better coupling between 
the source in front of the lens and especially for the most sub- 
wavelength modes. The resonant nature of the modes matches 
the impedance of the small electric dipole, or equivalently, the 
Purcell effect compensates for the weak radiation of the deep 
sub- wavelength modes. As a summary, we have seen that the 
Purcell effect implies a rate of spontaneous emission propor- 
tional to the quantity {Dk±)'^, while the radiation efficiency 
of a given mode is proportional to the inverse of this quantity. 
Overall, the two effects counterbalance themselves ensuring 
that all of the subwavelength details radiate the same amount 
of energy. And the higher the Bloch wavenumbers, the thinner 
the bandwidth carrying the energy. 

On a temporal point of view, the Q factor is the opposite 
of the attenuation time of the system: a higher quality factor 
implies a lower damping, and so high Q modes oscillate for 
longer times. This is why in the original paper of metalense^ 
we presented a curve of the lifetimes of the modes inside the 
structure versus their Bloch wavenumbers kj_. Here, we have 
explained why the data were growing with k± . The linear fit 
initially obtained results from the short time accessible for the 
time domain simulation compared to the extremely high Q 
factor attainable. 

In order to seal the validity of this physical description, we 
measured the return loss of a small dipole with a network an- 
alyzer (Agilent N5230C) with and without the presence of the 
lens (Fig. [6]). The return loss is the ratio of the power of the 
outgoing signal to the power of the signal reflected back from 
the small dipole. If the return loss is equal to 1, the whole 
input power is reflected back meaning that there is no radi- 
ated field. While when this quantity decreases, some amount 
of energy has been transmitted by the dipole: this energy 



10 



corresponds to the radiated energy. Without the lens, it re- 
sults in a flat curve in ampHtude demonstrating that the smaU 
dipole is totally reactive, that is to say, it does not radiate any 
field. In the presence of the resonant metalens, the curve de- 
creases down to a value of 0.7, meaning that the source is 
better matched due to the Purcell effect as stated before. The 
better coupling between the source and the resonant medium 
permits to compensate the less efficient radiation of the most 
sub wavelength modes: the two effects counterbalance them- 
selves and each mode eventually radiate the same amount of 
energy. 



keep the formalism of their article, it is convenient to use 
the notation tzju = ^Ttm- Also, by analogy, we denote the 
quantity k, = — k±'^ = ijo, where 70 is a real quan- 
tity for the waves we consider here (evanescent waves in free 
space). Then, an Additional Boundary Condition is intro- 
duced to solve the boundary problem at the air/wire medium 
interface^^, and the law obtained can be summarized in: 



(A.l) 



VI. CONCLUSION 



with. 



In this article, we have explained the physical mechanism 
of our original paper concerning the concept of resonant 
metalens ^ The intuitive property of resonance splitting when 
considering coupled resonators has been investigated through 
a wave approach. Starting from the wavefield inside an infi- 
nite wire medium, we have seen that a wire medium with fi- 
nite dimensions can be seen as a macro resonator that contains 
many resonance frequencies. Interestingly, each of these reso- 
nances corresponds to an eigen wavefield that fluctuates on its 
own sub- wavelength scale. Then, we have shown how these 
modes can radiate into the far field region information that is 
usually contained in the evanescent spectrum. Eventually, by 
placing such a resonant medium near a sub- wavelength source 
we have seen that the Purcell effect permits an efficient radia- 
tion. 

These physical mechanisms permit to get information of 
a source on a smaller scale than the wavelength, stored as a 
time/frequency signature, in the far field region. In the next 
article we will explain how one can use these properties for 
breaking the diffraction barrier in a focusing/imaging scheme 
and study the main limitations as well as losses impact. 



Appendix: Sophisticated dispersion law 

When calculating the dispersion law of guided modes in- 
side the wire medium, we have took sides for neglecting both 
the TM evanescent waves and the high orders of diffraction 
generated at the reflection. For the range of frequency con- 
sidered, the TM waves are strongly evanescent because the 
plasma wave number kp is really higher than any other wave 
number considered^ ^ . And the high order of diffractions are 
also strongly evanescent because the fundamental one is al- 
ready evanescent. All of these waves should modify the re- 
flection coefficient of equation ([12]) by an additional factor 
which affects phase, and thus dispersion. 

A dispersion law of the guided modes considering the TM 
waves is presented in the appendix of ref.^. To obtain this dis- 
persion law the authors used an effective medium approach 
which leads to an effective anisotropic permittivity along the 
wire axis^. It permits to obtain the Bloch wave number of 
the TM modes: hiju = \J^^ — — k^^, which is a com- 
plex quantity for the range of frequency considered here. To 



k_\_ -\- kri 



!4tanh(:^) (A.2 



7app = 7o ^ , '2 " + 7tm "tt ^^^^^ ) '^^■2) 



This dispersion relation is actually the same as the one pre- 
sented in the manuscript with a correction due to the TM 
waves which appears in 7app. In the limit of thin wires and 
for frequencies well below the cutoff of TM waves (which is 
the case considered here), the plasma wave number kp is ex- 
tremely higher than any other wave number considered. In 
that limit case we find that 7app = 70 and the dispersion law 
presented in equation ( 15 ) remains true. With the geometrical 
parameters we took in this example, the difference between 
the two laws is less than 1 % (Fig. [7]): this is the reason why 
we chose not to take into account the strongly evanescent TM 
waves. 



As seen in the expression of 7app which is always greater 
than 70, the TM waves tend to increase the apparent wave 
number associated to the evanescent decrease in air, or equiva- 
lently to decrease the penetration depth inside air. The guided 
modes are much more weakly confined to the interface than 
the case without considering TM waves. As a consequence, 
for a given transverse wavenumber k± it results in a higher 
frequency for the guided mode than the dispersion law pre- 
sented in the paper. 

We would like to argue that the dispersion law obtained 
with 3D time domain simulations (Fig. [7]) gives a distinct 
comportment for large k± than the one expected when adding 
the TM waves. Indeed, while the correction due to the TM 
waves tend to increase the frequency, the simulation results 
show that a decrease is expected. Actually, the periodic na- 
ture of the wavefield and thus the Bragg waves have been 
neglected, a remaining problem with effective medium ap- 
proach. Here we propose to do the calculation of the disper- 
sion law taking into account so me h igh orders of diffraction 
with a method inspired from ref 

We consider a TEM Bloch wave with transverse wave num- 
ber = /cxGx traveling in the half space z < filled out 
with a wire medium. This Bloch wave encounters a boundary 
with vacuum at z = and gives rise to a reflected TEM Bloch 
wave inside the wire medium region. The superposition of the 
2 TEM Bloch waves writes at z = 0: 
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To eliminate the unknowns, we use the fact that the tan- 
gential components of the electric field must be continuous 
at the vacuum/wire medium interface over the entire unit cell 
(which represents the area ft), while the magnetic field com- 
ponents are continuous only at the air regions of the unit cell 
(the area Q minus the area V corresponding to the wire sec- 
tion). Matching the x component of the electric field, multi- 

plying by e~*^^ ^, and integrating over the area \ t yields: 
where, 



FIG. 7. The dispersion laws of the first resonance frequencies 
fm (normalized to the original Fabry-Perot frequency /o) versus 
the transverse Bloch wave number k± (normalized to the original 
Fabry-Perot wave number ko) from equation ( p3] l only considering 
TEM waves in the wire medium (blue solid line), from ref ^ with 
our parameters considering TM and TEM waves (red dashed line), 
and from equations and \A.9\ considering only TEM waves 

but high diffraction orders (green dashed dotted line). Superimposed 
with stars is the results obtained from 3D simulations. 



(m) 



(A.6) 



Again, matching the y component of the magnetic field, 
multiplying by F(x, y)e~^^^^, and integrating over the area 
Q — V gives: 



t(-)/,* = (l-r)Jo 



(A.7) 



F{x,y) 

ETEM = (l + r)e^^-- ( 





1 







Htem = ^(1 - r)e''^^ F{x,y) \ (A3) 

MOC \ Q 

The function F{x^y) stands for the Bloch envelope of the 
TEM wave inside the wire medium, and notably this function 
is null over the wire section of a unit cell. One can notice that 
the y (respectively x) component of the electric (respectively 
magnetic) component of the TEM wavefield has been ne- 
glected: the diffraction orders along the ?/-axis are neglected. 
This implies that the reflected and transmitted TE modes will 
be neglected too, but in the general case they must be taken 
into account. We will also neglect the TM reflected waves 
because their influence has been discussed earlier. Thus, the 
transmitted field in the vacuum region, with these approxima- 
tions, writes at z = 0: 



Ev; 



Hv 



V -k'f^/k , 

- 2^ — e I 1 I (A-4) 



C ^ flQC \ Q 



and in the range of frequency consid- 



with k^^^ =k±^^ 
ered, we have again /ci^^ = ijn = ^j^^ — 



where, 



Jo 



F{x,y)HS 



(A.8) 



Then we can extract the reflection coefficient of the TEM 
wave, and as done in the manuscript when adding a second 
along the wires, it eventually gives the same dispersion law as 
equation \A.l\ with in this case: 



7app 



Jo 



E 



(A.9) 



To evaluate this quantity, we need to give an expression for 
the function F{x^y). For simplicity we have chosen that this 
function is equal to 1 in the domain uj — V, and null in the do- 
main D. The dispersion law obtained considering the 2 first 
diffraction orders (+1 and —1) is presented in figure [7] One 
can notice that this dispersion law gives a better fit to the re- 
sults obtained from the full 3D simulations. Obviously, this 
calculation is not an exact one but it has the merit to show 
that taking into account the diffraction orders tends to increase 
the apparent penetration depth inside air, thus giving a better 
comportment for the high transverse wave numbers. Never- 
theless, this calculation greatly complexifies the understand- 
ing for giving a correction less than 1 %: this is reason why 
we also neglect the high diffraction orders during the whole 
manuscript. 
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